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This letter considers the dynamics of a stiff filament, in particular 
the coupling of twist and bend via writhe. The time dependence 
of the writhe of a filament is ~ Lt^^^ for a linear filament 
and ~ t^l"^ /L for a curved filament. Simulations are used to 
study the relative importance of crankshaft motion and tube like 
motion in twist dynamics. Puller's theorem, and its relation with 
the Berry phase, is reconsidered for open filaments. 

DNA and other stiff polymer systems such as actin filaments are of interest both for their 
intrinsic biological importance, but also due to the fact that micro-manipulation techniques 
give a very detailed vision of their properties, allowing one to study fundamental processes 
in polymer statics and dynamics. Recent theoretical and experimental work on DNA force- 
extension curves has shown the importance of the static coupling between bend and twist 
in stiff polymers Non-trivial twist- writhe correlations in flexible polymers have been 

recently found in a simple dynamic lattice model in which linking number is conserved 0. 
Here I examine the dynamics of twist and writhe fluctuations in order to clarify the driving 
forces and dissipative processes important in the motion of filaments. I use a combination 
of scaling arguments and simulations; the numerical studies are complementary to recent 
work which considers the zero temperature motion of stiff filaments via integration of 
the continuum differential equations. Here I study Brownian dynamics of stiff filaments and 
introduce a generalized bead-spring model suitable for studying twist and writhe dynamics 
of polymers. 

There are two physical processes which can lead to the rotation of the end of a filament 
about its local tangent. The two mechanisms of spinning are twist and writhe. Twist 
corresponds to excitation of internal torsional degrees of freedom of the filament. Writhing 
is due to the three dimensional geometry of bending of the filament in space. We shall 
reserve the word "spinning" for this process of end rotation and denote the spinning angle 
by The two contributions are additive so that we can write \1/ = + Wr, with (f) the 
contribution from internal twisting and W^. the contribution from the geometry of the path 
in space. 

Excitations of internal twist degrees of freedom are well understood and have been 
observed in DNA and actin via depolarized light scattering experiments. The writhing of 
a filament is best expressed with the help of Fuller's theorem, [||. Consider a filament 
(with boundary conditions such that the two end tangents are maintained parallel) which 
is bent into a non planar curve. The tangent to the filament t(s) sweeps out a curve on the 
unit sphere, S2 as a function of the internal coordinate s. It can be shown, |]^, that the 
contribution to the spinning due to the writhe is just the area, Aq, enclosed by t(s) on ^2. 

Let us assume that we are working with uniform filaments with circular cross sections 
then the two process are controlled by two independent elastic constants, the torsional 
stiffness, K and the bending stiffness If we work in units such that ksT = 1 these 
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two elastic constants have the dimensions of lengths, and are indeed the persistence lengths 
for static torsional and bending correlations of the filament. The dynamics obeyed by the 
torsional and bending modes are, however, different which leads to distinct contributions to 
as we shall now demonstrate. 

The torsional motion of a straight filament obeys a Langevin equation 

where the angle t) is the local rotation of the filament in the laboratory frame, is 
the thermal noise, a is comparable to the radius of the filament, (we neglect all pre-factors 
of order unity and use units in which 77 = ksT = 1, in this system of units time has the 
dimensions of a volume, so that Is = 4/im^). The solution of eq. (|T]) is well known: It is 
characterized by dynamic correlations of the angle which diffuse along the filament according 
to the law l^^ ~ tK/a?. A monomer on the filament rotates by an angle (0^) ~ ^Jt/ Ko?^ 
if Itw < L. When Itw > L the twisting motion is sensitive to the boundary conditions on 
the filament; for free boundary conditions the filament rotates freely and {(p'jree) ~ t/La?. 
If one end is held at a fixed angle the other has an amplitude of rotation which saturates to 
(0Ld) - L/K. 

Let us now calculate the writhing contributions to the spinning. For a short filament 
{L/k < 1) oriented along the z-axis the area enclosed by the curve t(s) is given by = 
1/2 / 62.(1 A t)ds. The writhe of a filament is signed and averages to zero for a polymer at 
thermal equilibrium. We shall consider the statistics and dynamics of {A^[t)AQ{t')) . The 
linearized bending modes of the filament obey a Langevin equation 

^^-.^.^.(M) (2) 

with the transverse fluctuations of the filaments and the Brownian noise. Using eq. 
(0) we find that (tj, (g, 0)tj(— g, t)) = 6ij exp{—KqH) / nq'^ with {i,j} = {x,y}. The fourth 
order correlation function in {A^{t)AQ{t')) can be expanded using Wick's theorem: A short 
calculation shows that 

{AD = {W^) - L'/n\ (3) 

and that 

(^{Wr{t) - WM?) ~ hm/K" ~ t'l\ (4) 

with liit) ~ {ntY^^ the characteristic length scale in solutions of eq.(^. We note that the 
same result can be found in a scaling approach using the following argument: The path 
t(s) on ^2 is a Gaussian random walk ; On a time scale t each section of length li of the 
filament re-equilibrates a random walk with radius of gyration on ^2 of ~ and area 
Ai ~ ±/i/k. There are A'^i = L/li dynamically independent contributions to the writhe 
variation giving {/SW^) ~ NiA\ = li{t)L/n'^ as found above. This argument shows in 
addition that the scaling law eq. (^ is valid at short times even for filaments for which 
L> n. 
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What influence does this writhe fluctuation have on the global motion of a fllament. 
To answer this question consider two extreme, non-physical, cases before coming back to 
the typical experimental situation. Consider, flrstly, the case K/k ^ 1; spinning of the 
end of a polymer can only occur via writhe, the twist degrees of freedom are frozen out. 
The derivation of eq. has neglected rotational friction: It has been derived using a 
description of only the transverse motions of the polymer. It can only be valid when the 
rotational friction of a fllament is so small that spinning driven by the writhe is able to 
relax without build up of torsional stress. Similar remarks have recently been made in the 
propagation of tensional stress in stiff polymers [[^]. In this case simple arguments allow one 
to deduce the validity of the approximation made by balancing the driven motion due to 
the transverse fluctuations against any additional sources of friction. We shall now apply 
the same argument to the torsional motion, before checking the results numerically. For 
the rotational friction to be negligible the free rotational diffusion of a section of fllament of 
length L must be faster faster than the driven writhing motion, implying (^j^ee) > {^W^) or 
t/La^ > li{t)L/ n"^. For fllaments which do not satisfy this criterion, i.e. when L > Iwri^) = 
ylfK/a?, the torsional stress has not had time to propagate between the two ends of the 
filament and the law does not apply. If we observe now the end of a filament for times 
shorter than the torsional equilibration time the spinning is due to an end section of length 
Iwr and we should substitute this effective length in eq. to calculate the end motion 
(AW^) ~ fi"^ /an^/'^ ~ t^/^. It is only after the propagation of the torsional fluctuations 
over the whole length of the filament that (D becomes true. We thus hypothesize a scaling 
behavior for the spinning of a filament with K/ k ^ 1 

^ ' ft'/' ^ l,o8/3L8/3 ) • ^ ' 

with Q{x) ~ x^^^ for x small and Q{x) ~ 1 for x large. 

A second extreme case is = 0, as is for instance the case for many simple numerical 
bead-models of worm-like chains. In this case there is no need for the beads to rotate to 
follow writhe, the writhe fiuctuations are absorbed by the internal twist degree of freedom 
without cost. The writhe does not lead to spinning; we do expect however that there are 
strong dynamic (anti) correlations between twist and writhe over the length of the filament. 
We conclude that the law is valid for all times, it has no consequence however on the 
rotational dynamics of a monomer on the chain that could be detected, for instance, in a 
scattering experiment. 

Let us now turn to the physical case, K/k ~ 1. From the expressions for and kw we 
find that Uw/lwr = ^ <^ 1 at short times. On the length scale Itw the twist field and the 
writhe of the filament are in equilibrium and we can add the fiuctuations of the twist and 
writhe degrees of freedom. Beyond Itw the successive section of the polymer are dynamically 
decoupled: a writhe fiuctuation can not be transmitted faster than the signal transmitted by 
Itw and the second scenario becomes vahd with Itw playing the role of a bead size. Beyond Itw 
the filament can writhe freely without any local consequence on the dynamics of the chain. 
Thus the experimental consequences of the writhe fiuctuations are probably negligible for 
times such that Itw < L. When Itw reaches L the amplitude of twist fiuctuations saturates, 
while the dynamics of the slower bending modes continues to evolve with varying as eq. 
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(I) until li{t) = L. These additional fluctuations could be observed experimentally. 

These arguments have been checked by simulations on a discretized bead-spring model: 
Take iV + 1 spherical beads of diameter a connected by stiff harmonic springs. Each bead 
is characterized by its position r and a triad of orthonormal vectors. M = {b, n, t}. The 
vector t is an approximation to the local tangent to the filament at r. {n, b} span the 
normal space of the filament; n could describe the position of some feature, such as the 
large groove in DNA on the surface of the molecule. The springs are not connected to the 
centers of the beads rather a bead i is linked to its neighbors by connections situated on its 
surface at {r^ — atj/2, + atj/2}. The energy is 

E = B/2 J2{r^ - r.+i - ati/2 + aU+i/2f 

i 

i 

+ ^^740 ^(1 - n, ■ rii+i - h, ■ bi+i + U ■ U+i) (6) 

i 

The first term imposes the linear topology of the filament; when B is large the curvilinear 
length of the filament is L = Na. In the simulations we are only interested in the limit of 
large B. The second term is a bending energy. The third term is the torsional energy. I 
integrate the equations of motion for the position and orientation adding friction and thermal 
noise coupled to the translational and angular velocities, effectively leading to Brownian 
dynamics. Since we are particularly interested in the writhe dynamics I use boundary 
conditions where the positions of the ends of the filament are free to move but the end 
tangents are constrained via external torques, so that the writhe can be calculated via 
Fuller's result. 

A first series of simulations were performed to study the internal twisting dynamics of 
very long filaments. This simulation was performed to confirm the intuition that twist 

modes can transport torsional stress even in the presence of large bends in the filament. This 
stress is carried by rapid spinning of the filament in a slowly evolving tube, rather than the 
rotation of the whole filament collectively moving against solvent friction, (like the rotation 
of a rigid crankshaft). Simulations are performed for filaments of varying length with two 
possible boundary conditions. The first boundary condition is that both ends are free to 
spin and we look for a crossover from the dynamics of internal modes described by eq. 
(|^) to free rotation. In the second case we block the rotation of one end and check that 
this leads to a saturation of rotation angle. The scahng displayed by Fig. (P shows that 
torsion propagates perfectly well over many persistence lengths without hindrance from the 
tortuosity of the path in space. Eq. (|lD is valid for thermally bent filaments as well as for 
torsional fluctuations around a linear configuration. 

A second family of simulations were performed to study the writhe and twist fluctuations 
of filaments both with very high torsional constants to verify the predictions made for the 
writhe dynamics. Fig(2) shows the the scaling proposed above, eq (|^), is indeed seen in the 
case of high torsional stiffnesses. 

In order to calculate the writhe of a curve via Fuller's theorem we have until now consid- 
ered paths t(s) which are closed on 1S2. In general this is most inconvenient: there are many 
experimental situations where one would like to compare the spinning of filaments oriented 
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FIG. 1. Spinning of an end as a function of time. Two different boundary conditions are applied 
at the opposite end: free spinning or held. At short times the boundary conditions are unimportant. 
The angle evolves as (i;^^) ~ t^/"^. For long times there is either saturation of the twist angle to 
i^hoid) (lower curves) or free diffusional rotation of the whole filament according to {(p'^j.^f,) (upper 
curves), k = 30a, L varies from k/3 to 6.6k. K/k = 1. Splitting of lower curves is due to writhe. 



in an arbitrary direction. We now consider how one might generalize the idea of writhe and 
spinning to filaments with arbitrary open boundary conditions. In general it is impossible 
to compare the two rotation frames M(L) and M(0). The main ambiguity comes from the 
fact that rotations are non-commutative. However for filaments for which L/k, the non- 
commutative nature of the rotations is a higher order correction and we can consider that 
in going along the filament oriented along z we have a rotation which can be unambiguously 
decomposed into its three Cartesian coordinates {fix, fiy,fiz)- We shall proceed by studying 
equations d'M{s)/ds = 0(s)M(s) where Oi j = eijkUJ^ with oj'^ = dtk/ds for k = {x,y}. 
The angular velocities uj^ and a;^ correspond to bending the filament, = in the absence 
of twist. We iteratively integrate this equation. The lowest order contribution to rotations 
about the z-axis comes in second order. 

n, = l f' dsds' d{s - S') {uj:ujI - U;:,UJI). (7) 

2 JO 

Integrating by parts to transform the 6'-function into a (5-function gives 

Vt, = ]^e,.l^j\iM)ds + t{Q) M{L)^ (8) 

For a closed curve on ^2 the first term is the area enclosed by the curve t(s) on a patch of 
a sphere in agreement with Fuller's theorem. The second, boundary, term present when the 
curve is open corresponds to closing the path by a geodesic from t(L) to t(0). Thus we have 
a simple generalization for Fuller's theorem valid for short chains allowing one to calculate 
the writhe state of a polymer with fluctuating boundary conditions. 
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FIG. 2. Left: Spinning of the end of a linear filament with K/ k = 300. Data is scaled according 
to eq. (5). Right: Spinning of the end of a semicircular section of filament for two values of K. In 
the lower curve all spinning is due to writhe. For the upper curve there is a two time behavior, 
with firstly twist then writhe relaxing. Simulations are for L/a = k/o = 35. 



We remark that Fuller's theorem is closely related to Berry's phase in experiments on 
polarized light transmission along bent optical fibers |]T2[. The vector t(s) corresponds to the 
local tangent to the fiber, while the vectors, {n, b} transform in exactly the same manner 
as the plane of polarization, via parallel transport on 1S2. The problem of closing paths 
in 1S2 has close analogies in the treatment of Berry's phase in non cyclic Hamiltonians. In 
quantum and optical systems interference phenomena allow one to define the relative phase 



of a system even if the evolution has not been cyclic |0 via the Pancharatnam connection. 
This convention also closes an open trajectory again via a geodesic; in wave physics this 
convention is even generally applicable and is not limited to small paths in the appropriate 
projective Hilbert space. 

All the above discussion has been for a filament for which the round state is linear. Let us 
repeat the arguments leading to (|) for short bent filaments to study the writhe fluctuations 
of short curved DNA sections. The pre-existing bend considerably modifies the arguments. 
Consider a section of filament which bends an angle $ at zero temperature due to intrinsic 
curvature. At zero temperature the bent loop has a tangent map which maps to the equator 
of ^2. Under the map to ^2 a distance of s in real space becomes ^s/L. At finite temperature 
there is competition between this constant drift and the thermal agitation which leads to a 
Gaussian random walk on ^2. Let us now consider the dynamics in a scaling picture: For 
the very shortest times only small scale structure is changing:- the writhe dynamics reduces 
to the case discussed above. At longer times the stretched path moves collectively:- each 



section of length 5i = on ^2 moves up or down by = yli/n. Counting Ni = L/li 

dynamically independent sections 

/\W^ ~ Ni{5irnf = ^HI/kL ~ t^/^. (9) 
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The exponent for the writhe fluctuations has changed from 1/4 to 1/2 due to the stretching 
of the configuration in 1S2. This prediction is tested in fig. (2), where we have added a 
spontaneous bending energy, in aSjiii.ti+i to curve the filament. We now see two distinct 
regimes in \E'^ ~ t^/^ due to the intrinsic twist dynamics and then the writhe dynamics, the 
writhing contribution is also much enhanced over its value for linear filaments. It would also 
be interesting to study the case of closed loops, but the simple scaling argument used here 
becomes considerably more complicated due to the non-local integral constraint J t{s) ds = 
coming from closure. 

This paper has studied the short time regime dominated by the bend and twist rigidities, 
rather than the very long time dynamics where crossovers to Zimm or Rouse dynamics 
occur. Our arguments have been for the simplest model of a stiff polymer, of uniform 
cross section and without disorder. The Marko-Siggia energy function []T^ already has rich 
static behavior and one would anticipate that the cross terms in this Hamiltonian could 
increase the importance of the dynamically effects discussed here. Similarly disorder in the 
ground state is expected to dramatically modify the dynamics, by giving a preformed writhe 
on ^2 and by modifying the picture of simple spinning of the polymer in its tube due to 
new dissipative processes; It has recently been argued that a mixture of crankshaft and tube 
spinning should co-exist even in the case of weak disorder coming from sequence fluctuations 
in DNA 0. 

I would like to thank A. Ajdari, R. Everaers, F.Julicher, P. Olmstead, C. Wiggins and 
T.A. Witten for discussions. e 



[1] C. Bouchiat and M. Mezard. Phys. Rev. Lett. 80, 1556, (1998). J. D. Moroz, P. Nelson. 

Macromol. 31, 6333, (1998). 
[2] T.R. Strick et al. Science 271, 1835 (1996). 
[3] J.S. Plewa, T.A. Witten. Preprint 

|cond-mat/9909367 

[4] R. E. Goldstein, T.R. Powers, C.H. Wiggins. Phys. Rev. Lett, 80 5232-5235, (1998). 
[5] M.D. Barkley, B.H. Zimm. J. Chem. Phys., 70:2991-3007, 1979. 
[6] F.B. Fuller. Proc. Natl. Acad. Sci. 68, 815, (1971), 75, 3557, (1975). 
[7] R. Everaers et al. Phys. Rev. Let. 82, 3717-3720, (1999). 
[8] R.D. Kamien. Eur. Phys. J. B, 1 1, (1998). 
[9] P. Nelson preprint. |cond-mat/9907040| 
[10] J.M. Schurr et al. J. Chem. Phys 106 815 (1997). 

[11] C. Levinthal, H.R. Crane. Proc. Natl. Acad. Sci. USA, 42 436-8, (1956). 

[12] F.D.M. Haldane. Opt. Lett. 11 730, (1986). 

[13] J. Samuel, R. Bhandari. Phys. Rev. Lett. 60 2339, (1988). 

[14] J.F. Marko, E.D. Siggia. Macromolecules 27, 981, (1994). 



7 



